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We  show  that,  in  addition  to  the  known  monohydrate,  LiOH  forms  a  dihydrate  at  elevated  pressure. 
The  dihydrate  involves  a  large  number  of  H-bonds  establishing  chains  along  the  (001)  direction. 
In  addition,  the  energy  surface  exhibits  a  saddle  point  for  proton  locations  along  certain  O  inter¬ 
atomic  distances,  a  feature  characteristic  for  superprotonic  conductors.  However,  MD  simulations 
indicate  that  LiOH  •  2H2O  is  not  a  superprotonic  conductor  and  suggest  the  relevant  interpolyhe- 
dral  0-0  distances  being  too  large  to  allow  for  proton  transfer  between  neighboring  Li-coordinated 
polyhedra  at  least  on  the  time  scale  of  the  MD- simulations.  ©  2011  American  Institute  of  Physics. 
[doi:  10. 1063/1. 3543797] 


I.  INTRODUCTION 

Liquid  H2O  is  a  solvent  for  many  polar  compounds; 
however,  at  ambient  pressure,  water  ice  accepts  only  trace 
amounts  of  chemical  impurities.^  However,  at  elevated  pres¬ 
sure  or  at  low  temperatures,  molecules  or  atoms  of  van  der 
Waals  radii  much  larger  than  that  of  H2O  form  hydrous 
compounds,  where  hydrogen-bonded  ice  networks  encapsu¬ 
late  these  molecules  in  large  cages  of  so-called  clathrate 
structures'’^  including  clathrates  of  H2."^“^ 

Recently,  it  was  found  that  ice  VII  dissolves  substantial 
amounts  of  LiCl.^  These  hndings  inspired  us  to  search  for 
light  element  hydrides  forming  clathrates  or  zeolitelike  com¬ 
pounds  with  water  ice. 

Under  ambient  conditions,  LiH  reacts  with  water  to  form 
free,  gaseous  H2  and  LiOH.  LiOH  by  itself  can  form  a  crys¬ 
talline  ionic  monohydrate. ^  While  we  observe  that  similar  re¬ 
actions  prevail  at  higher  pressures,  at  least  up  to  3.5  GPa,  we 
hnd  that  LiOH  forms  a  dihydrate  at  elevated  pressure,  which 
exhibits  features  characteristic  for  proton  conductors:  Li  co¬ 
ordination  polyhedra  that  can  easily  rotate  and  undersaturate 
coordinating  oxygen  atoms  that  can  bind  hydrogen. 

II.  EXPERIMENTAL 

Single  hexahedral  crystals  of  LiH  (Alfa)  with  50  ptvc? 
edge  lengths  were  loaded  individually  in  the  sample  chambers 
of  a  diamond  anvil  cell.  The  chamber  size  was  chosen  such 
that  half  of  the  chamber  was  occupied  by  the  hydride  crystal. 
The  chamber  wall  was  established  by  a  Re  gasket  of  100  /xm 
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thickness  and  with  central  bore  of  150  /xm  diameter  and  the 
culet  of  one  diamond  anvil.  A  spherule  of  annealed  ruby  of 
5  /xm  diameter  was  placed  close  to  the  inner  border  of  the 
chamber.  The  second  diamond  anvil  was  then  placed  near 
the  top  of  the  gasket  surface,  leaving  a  gap  of  several  10  /xm. 
Deionized  water  was  injected  with  a  syringe  into  this  gap 
between  anvil  and  gasket,  while  the  whole  assembly  was  ob¬ 
served  under  a  microscope.  Once  water  hhed  the  chamber,  the 
upper  anvil  was  lowered  within  1  s  to  close  the  chamber  and 
to  establish  an  initial  pressure,  which  was  clearly  beyond  the 
fluid  H2-water  coexistence  boundary  at  300  K  since  no  gas 
formation  was  observed  at  this  initial  pressure.  The  mechani¬ 
cal  load  was  raised  subsequently  to  the  experimental  pressure. 
Three  loadings  were  prepared  at  1.5,  2.0,  and  3.5  GPa. 

Ah  samples  were  examined  by  x-ray  diffraction  at 
the  high  flux  microdiffraction  undulator  beamline  ID-B  at 
sector  16  (HPCAT)  of  the  Advanced  Photon  Source  (APS). 
A  34.325  keV  x-ray  beam  was  focused  in  vertical  and 
horizontal  directions  by  Kirkpatrick-Baesz  mirrors^  to  a 
5x5  /xm^  spot  at  the  sample  chamber.  The  sample  was 
rastered  in  a  plane  perpendicular  to  the  beam  direction  over 
a  100  X  100  /xm^  area  during  signal  accumulation  in  order 
to  maximize  crystallite  powder  statistics.  Diffracted  radiation 
was  recorded  with  an  MAR345  image  plate  detector  for 
average  accumulation  times  of  300  s  per  frame.  The  recorded 
patterns  were  corrected  for  geometric  distortion  and  in¬ 
tegrated  along  the  azimuthal  angle.  Instrumental  profile 
parameters  were  extracted  from  Le  Bail  refinements^ of 
the  patterns  of  an  NIST  Ce02  standard. 

In  all  experiments  we  identified  LiOH  •  H2O,  ice  VI,  ice 
VII  (at  2.0  and  3.5  GPa),  and  hydrogen-hydrates  Cl  (at  the 
lowest  pressure)  and  C2.^  The  presence  of  lithium  hydrox¬ 
ide  monohydrate  and  of  hydrogen-clathrates  in  ah  three  runs 
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FIG.  1.  Calculated  (line)  and  observed  (crosses)  powder  diffraction  pat¬ 
tern  of  LiOH  •  2H2O,  phase  C2,  ice  VII  and  LiOH  •  H2O  at  3.5  GPa.  Tick- 
marks  indicate  the  2?^  angles  of  allowed  reflections  of  LiOH  •  2H2O,  C2,  and 
LiOH  •  H2O,  respectively.  The  bottom  line  represents  the  residual  of  calcu¬ 
lated  and  observed  pattern.  All  signiflcant  deviations  (e.g.  around  15°  2i^) 
are  from  misfit  of  the  C2-phase  (see  experimental  section  for  details). 


indicates  that  LiH  had  reacted  with  water.  Free  hydrogen  was 
captured  in  the  Cl  and  C2  clathrates.  No  LiH  was  identified 
in  any  pattern. 

In  the  run  at  2.0  GPa,  we  identified  ice  VII,  hydrogen- 
hydrate  Cl,  and  LiOH  •  H2O.  In  the  run  at  3.5  GPa,  we  found 
hydrogen-hydrate  C2,  ice  VII,  LiOH  •  H2O,  and  the  diffrac¬ 
tion  pattern  of  an  unknown  phase.  The  pattern  does  not  belong 
to  any  known  phase  of  ice,  Li-oxide,  -hydroxide,  or  -hydride. 
The  most  intense  diffraction  features  occur  between  7.7°  and 
10.3°  2z^,  while  there  is  a  high  density  of  reflections  of  low 
to  moderate  intensity  at  higher  angles  (Fig.  1)  with  increasing 
overlap  at  higher  angles.  As  we  will  show,  this  new  phase  is 
LiOH .  2H2O. 

The  remaining  large  number  of  reflections  indicates  pres¬ 
ence  of  a  rather  low  symmetric  structure  or,  alternatively,  a 
mixture  of  unidentified  phases.  Besides  the  intrinsic  difficulty 
of  determining  an  unknown  structure  from  powder  data  col¬ 
lected  from  a  sample  confined  in  a  diamond  anvil  cell,  we  face 
the  additional  problem  of  unknown  stoichiometry.  We  started 
to  examine  this  unidentified  part  of  the  diffraction  pattern  by 
comparison  to  calculated  patterns  of  known  ice,  clathrate,  and 
zeolithe  structures.  We  also  calculated  Li-bearing  ice-network 
structures  using  known  Li-silicate  structures  and  substitut¬ 
ing  Si  and  A1  by  O  and  O  by  protons  and  subsequent  com¬ 
pression  of  those  model  structures  to  densities  comparable 
to  LiOH  •  H2O  or  ice  VII  within  30%  margins.  This  search 
showed  similarity  between  the  observed  pattern  and  a  calcu¬ 
lated  ice  network  derived  from  spodumene.^"^ Another  approx¬ 
imate  match  was  found  for  a  layered  structure  of  Li  alternating 
with  a  two-dimensional  ice  network,  which  to  our  best  knowl¬ 
edge  has  no  immediate  analog  among  known  structures.  We 
used  these  initial  guesses  as  starting  sets  for  reversed  Monte 
Carlo  (rMC)  modeling  using  Free  Objects  for  Crystallography 
(FOX).^^  All  rMC  runs  converged  to  one  structural  arrange¬ 
ment  of  atoms.  Several  indications  support  the  conclusion  that 
the  result  of  the  rMC  modeling  is  not  arbitrary:  (a)  the  struc¬ 


ture  has  reasonable  interatomic  distances,  (b)  the  same  struc¬ 
ture  is  found  independently  on  the  starting  arrangement  of 
atoms  and  even  after  reducing  the  initial  set  of  refiections, 
(c)  although  the  rMC  calculations  were  run  without  symme¬ 
try  constraints,  the  resulting  structure  has  intrinsic  symmetry 
(see  below),  (d)  no  bias  by  interatomic  potentials  had  to  be 
used,  the  resulting  structure  is  robust  with  respect  to  the  ob¬ 
served  refiections,  and  finally,  (e)  additional  rMC  calculations 
using  Endeavor  1.1  (Ref.  16)  gave  the  same  result. 

An  automatic  symmetry  search^^  yielded  monoclinic 
cells  of  space  groups  Clcl  and  C12i/mL  However,  the  calcu¬ 
lated  pattern  for  the  Clcl  cell  matched  the  observed  pattern 
better. 

III.  COMPUTATIONAL  METHODOLOGY 

We  augmented  our  experiments  with  first-principle  den¬ 
sity  functional  theory  (DFT)  calculations  with  two  goals  in 
mind:  (1)  investigate  the  stability  of  the  structure  obtained 
from  x-ray  diffraction  and  (2)  to  perform  ab  initio  molecu¬ 
lar  dynamics  to  investigate  the  likelihood  of  hydrogen  mo¬ 
bility  below  the  melting  temperature  of  LiOH  •  2H2O.  All 
calculations  were  performed  by  using  the  software  pack¬ 
age  VASP.^^“^^  An  energy  cutoff  of  600  eV  was  found 
sufficient  to  obtain  converged  results.  The  electronic  band 
structure  was  determined  using  the  tetrahedron  method  with 
Bloechl  corrections.^^  The  interactions  between  ions  and  elec¬ 
trons  were  described  within  the  projector- augmented- wave 
formalism;^^’^"^  the  reference  configuration  (core  radii,  in 
units  of  Bohr  radii)  were  ls^2s^  (2.050  ae),  2s^2p'^  (1.520 
ae),  and  Is^  (1.100  ae),  for  Li,  O,  and  H,  respectively.  Elec¬ 
tronic  exchange  and  correlation  effects  were  treated  at  the 
GGA  level  in  the  parametrization  of  Perdew  et  During 
the  relaxation,  all  degrees  of  freedom,  lattice  and  atomic  po¬ 
sitions,  were  relaxed  simultaneously. 

Ab  initio  MD  simulations  were  performed  for  a  36  atom 
simulation  cell  of  composition,  4LiOH  •  2H2O.  The  cell  shape 
of  the  DFT  optimized  static  lattice.  The  temperature  was  con¬ 
trolled  by  a  single  Nose  thermostat.^^’^^  The  time  step  for 
the  integration  of  Newton’s  equation  of  motion  for  the  ions 
was  0.5  fs.  The  system  evolution  was  followed  up  to  4000 
timesteps  (2  ps).  We  analyzed  the  time  dependence  of  the 
mean- square  displacement  as  a  measure  of  mobility  for  that 
was  analyzed  for  Li,  O,  and  H  separately.  For  a  conventional 
superprotonic  conductor  it  would  be  expected  that  the  Li  re¬ 
mains  comparatively  stationary,  while  the  oxygen  and  the 
hydrogen  subsystems  should  show  larger  excursions  and  oc¬ 
casional  hydrogen  jumps  between  Li04  polyhedra.  Thus,  the 
mean-square  displacement  is  a  suitable  measure  for  the  as¬ 
sessment  of  mobility  in  superpro tonic  conductors.^^ 

We  compared  the  diffraction  pattern  of  the  relaxed  DFT 
structures  with  the  experimental  diffraction  pattern  (Fig.  1). 
The  experimental  pattern  exhibits  spotty  diffraction  features 
from  the  comparatively  coarse  grained,  coexisting  mono¬ 
hydrate.  These  nonstatistical  powder  diffraction  intensities 
affect  the  convergence  of  the  Rietveld  refinement.  Therefore, 
we  choose  a  different  approach  to  compare  experiment  and 
theory,  which  bases  on  clearly  identified  refiections  of  the  new 
phase  only  and  which  allows  for  removing  refiections  which 
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TABLE  I.  List  of  observed  and  calculated  structure  factor  moduli.  First  col¬ 
umn:  Miller  indices  of  the  observed  reflections,  second  column:  multiplicity, 
third  column:  observed  structure  factor  moduli,  and  fourth  column:  calculated 
structure  factor  moduli  based  on  the  ab  initio  calculated  model  structure.  The 
calculated  moduli  are  corrected  for  multiplicity,  while  the  observed  moduli 
represent  the  combined  values  from  contributions  of  overlapping  reflections. 


hkl 

Multiplicity 

IFclcI 

IFobsl 

1 

1 

0 

4 

4.33 

4.31 

1 

1 

-  1 

4 

1.82 

3.641 

2 

0 

0 

2 

28.88 

35.71 

0 

2 

0 

2 

5.67 

3.60 

1 

1 

1 

4 

46.14 

33.15 

0 

2 

1 

4 

4.67 

8.03 

2 

0 

-2 

2 

56.32 

62.96 

3 

1 

-  1 

4 

45.94 

45.65 

1 

1 

-2 

4 

28.80 

26.52 

2 

2 

0 

4 

37.38 

32.69 

0 

0 

2 

2 

3.05 

14.49 

3 

1 

0 

4 

15.10 

18.10 

3 

1 

-2 

4 

19.00 

21.70 

2 

2 

-2 

4 

37.38 

32.69 

1 

3 

-1 

4 

50.23 

47.71 

4 

0 

-2 

2 

4.58 

1.58 

0 

2 

2 

4 

28.68 

18.87 

2 

2 

1 

4 

7.41 

6.03 

1 

1 

2 

4 

14.29 

13.75 

1 

3 

1 

4 

5.2 

9.65 

3 

1 

-3 

4 

4.21 

3.96 

4 

2 

-2 

4 

18.97 

29.73 

5 

1 

-2 

4 

14.37 

14.44 

2 

4 

0 

4 

26.93 

22.10 

1 

3 

2 

4 

16.21 

12.97 

4 

2 

-3 

4 

27.56 

29.16 

2 

2 

2 

4 

5.62 

9.81 

5 

1 

-3 

4 

18.67 

24.37 

2 

4 

-2 

4 

5.85 

9.90 

are  affected  by  contributions  by  overlapping  Bragg  peaks 
of  coexisting  phases:  we  extracted  IF(hkl)rs  by  fitting  peak 
profiles  of  the  experimental  pattern  with  the  Le  Bail  method^^ 
and  compared  them  to  the  calculated  structure  factor  moduli. 
Hereby,  we  did  not  partition  IF(hkl)rs  of  overlapping  reflec¬ 
tions  but  compared  the  observed  IF(hkl)l’s  to  the  calculated 
ones  corrected  for  multiplicity.  Up  to  a  maximum  Q  of 
5.6  sample  reflections  could  be  distinguished  from  re¬ 
flections  of  the  coexisting  phases.  We  used  the  Wilson  plot  to 
identify  and  remove  IFobsl  strongly  affected  by  contributions 
from  coexisting  phases.  The  Wilson  plot  gives  an  average 
Biso  of  0.048(8).  Out  of  46  symmetry  allowed  reflections 
within  a  Q-range  up  to  5.6  there  were  29  reflections  of 
observable  intensity,  which  did  not  substantially  overlap  with 
diffraction  peaks  from  the  coexisting  phases  and  which  could 
be  used  for  evaluation  of  the  structure  model.  The  calculated 
and  observed  IF(hkl)rs  are  listed  in  Table  I. 

Based  on  this  set  of  reflections  we  calculate  a  relia¬ 
bility  factor  R=  E  (||Fobs -IFcidl/^ IFobsl )  of  0-17,  which  in¬ 
dicates  that  the  structure  calculated  by  DFT  represents  the 
structure  of  the  experimentally  observed  phase  at  least  to 
the  level  of  equal  bond  topology  and  coordination,  while 


fractional  coordinates  may  be  different  within  the  over¬ 
all  17%  margin  of  uncertainty  given  by  the  limited  qual¬ 
ity  of  the  experimental  data.  However,  comparison  of  the 
ah  initio  calculated  fractional  coordinates  with  those  obtained 
from  the  rMC  modeling  reveals  markedly  larger  deviations 
for  the  y-coordinates  than  for  the  x  and  z  coordinates.  This 
larger  discrepancy  could  reflect  an  unfitted  preferred  orienta¬ 
tion.  More  likely  it  is  the  result  of  omission  of  reflections  due 
to  overlap  with  the  diffraction  signal  of  the  coexisting  phases. 
In  sum,  the  Rint  of  17%  reflects  an  agreement  between  the  ex¬ 
perimental  and  ah  initio  calculated  fractional  coordinates  of 
3%-15%  for  the  x  and  z  coordinates  and  only  30%  for  the 
y-coordinate.  Since  the  ah  initio  calculation  converged  sta¬ 
bly  it  is  not  likely  that  the  y-coordinates  of  the  atomic  posi¬ 
tions  in  the  observed  material  deviate  much  more  from  the 
calculated  values  than  the  x  and  z  coordinates.  Fractional  co¬ 
ordinates  and  cell  parameters  are  given  in  Table  I,  and  the 
structure  model  is  shown  in  Figs.  2(a)  and  2(b)  in  different 
projections.  Upon  compression  or  decompression,  the  quality 
of  the  diffraction  patterns  of  lithium  hydroxide  dihydrate  de¬ 
grades  significantly  by  peak  broadening,  which  discouraged 
an  examination  of  the  structural  evolution  of  this  material  with 
pressure. 

IV.  DISCUSSION 

We  discuss  LiOH  •  2H2O  based  on  the  structure  calcu¬ 
lated  by  DFT.  All  atoms  reside  on  special  positions.  01  and 
02  have  two  protons  each  within  0.97-1.06  A  distance.  03 
has  only  one  proton  within  1 .02  A.  Li  has  the  closest  distance 
to  02  (1.94  A),  equal  to  the  Li-0  distance  in  LiOH-H20 
(Ref.  8),  while  all  other  0-atoms  are  not  closer  to  Li  than 
2  A.  Clearly,  01  and  02  belong  to  water  molecules,  while  03 
belongs  to  the  hydroxide  group. 

01  is  tetrahedrally  coordinated  by  H-bonds  to  adjacent 
01  and  02  atoms  (Fig.  1).  These  H-bonds  between  01  atoms 
form  a  chain  along  (100)  [Fig.  2(a)].  Another  set  of  0-H. . .  O 
interatomic  distances  involving  02  and  03  are  longer  than 
those  expected  for  H-bonds  (2. 9-3.0  A).  If  H-bonds  would 
form  along  these  directions  LiOH  •  2H2O  would  be  compati¬ 
ble  with  the  Pauling  rules  of  an  ice-network  forming  channels 
along  (010).  However,  the  static  position  of  the  proton 
attached  to  02  is  displaced  from  this  potential  bond  direction 
along  02-03,  making  the  0-H-H  distance  markedly  larger 
than  expected  for  H-bonds  and  distorting  the  hypothetical 
H-O-H  bond  angle  substantially.  This  displacement  is  likely 
due  to  the  influence  of  the  electric  field  of  Li,  which  prohibits 
the  formation  of  an  ice  network,  at  least  at  low  temperatures. 
The  Li-03  bond  vectors  and  the  03-H  bonds  establish 
together  a  fivefold  coordination  of  03,  while  Li  is  threefold 
pyramidal  coordinated  by  01,  02,  and  03  [Fig.  2(b)]  con¬ 
sidering  Li-0  distances  below  2.06  A  (with  a  fourth  distance 
at  2.33  A).  The  projection  in  Fig.  2(b)  shows  how  Li  resides 
between  the  H2O  tetrahedral  “chains.”  It  is  remarkable  that 
the  distances  between  Li  and  the  oxygen  of  the  hydroxide  on 
one  side  and  the  02  of  the  water  molecules  on  the  other  side 
differ  by  only  0.04  A,  making  Li  appearing  like  a  “network 
former.”  This  geometric  configuration  of  Li-0  and  H-0 
bonds  may  permit  for  efficient  proton  conduction  at  elevated 
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(b) 


FIG.  2.  (a)  Representations  of  the  structure  of  LiOH  •  2H2O  as  a  ball-stick 
model,  (a)  Projection  along  (100).  Blue:  O,  green:  Li,  white:  H.  White  lines 
indicate  the  directions  of  H-bonds.  The  most  remarkable  feature  of  the  struc¬ 
ture  is  the  arrangement  of  water  molecules  along  (100)  to  a  chain  of  four- 
member  rings.  The  lack  of  tetrahedral  coordination  of  all  oxygen  atoms 
but  01  violates  the  ice-rules  and  prevents  macroscopic  disorder  of  protons, 
(b)  Atom-  and  bond  colors  as  above.  Green  lines:  Li-O  distances  below 
2.0  A,  which  are  compatible  with  significant  ionic  interaction.  The  projec¬ 
tion  shows  how  Li  resides  between  the  H2O  tetrahedral  “chains.”  We  note 
that  the  distances  between  Li  and  the  oxygen  of  the  hydroxide  on  one  side 
and  the  02  of  the  water  molecules  are  shorter  by  only  0.04  A,  making  Li 
appearing  like  a  “network  former.”  However,  the  MD  simulations  indicate 
that  LiOH  •  2H2O  does  not  behave  like  an  itinerant  ice  network  and  does  also 
not  permit  itinerant  proton  diffusion  paths  which  would  make  this  material  a 
superb  proton  conductor. 

temperature,  if  correlated  Li-  and  0-motion  serves  as  a 
mediator  for  itinerant  proton  diffusion  paths. 

Elevated  temperature  in  combination  with  an  applied 
electric  potential  could  make  LiOH  •  2H2O  a  proton  conduc¬ 
tor  of  very  high  proton  density  but  rather  low  bulk  density. 


TABLE  11.  Atomic  positions  for  Li0H-2H20  calculated  ab  initio  at 
9.5  GPa  and  0  K.  The  calculated  cell  parameters  are  a  =  6.3714  A, 
b  =  7.6916  A,  c  =  5.0619  A,  and  ^  =  121.369°.  The  experimental  cell 
parameters  at  3.5  GPa  were  found  to  be  a  =  7.697(1)  A,  b  =  6.624(7)  A, 
c  =  4.518(4)  A,  and  ^  =  112.17(9)°.  The  space  group  is  Clcl  (#9). 


Wyckoff 

X 

y 

z 

Li 

4a 

0.3661 

0.1099 

0.3309 

01 

4a 

0.2754 

0.5574 

0.4952 

02 

4a 

0.0399 

0.2318 

0.1938 

03 

4a 

0.3026 

0.1246 

0.9144 

HI 

4a 

0.9421 

0.1247 

0.1055 

H2 

4a 

0.2837 

0.3069 

0.9673 

H3 

4a 

0.2467 

0.5370 

0.6673 

H4 

4a 

0.1303 

0.1112 

0.7500 

H5 

4a 

0.9531 

07059 

0.7908 

Based  on  the  mean- square  displacements,  the  molecu¬ 
lar  dynamic  calculation  at  1000  K  indicates  that  all  atomic 
species  move.  The  material  is  therefore  likely  molten.  At  500 
and  750  K  we  find  no  translational  motion  of  any  component 
and  atoms  move  about  their  equilibrium  positions.  In  partic¬ 
ular,  none  of  the  hydrogen  atoms  moves  within  a  2  ps  time 
period.  This  indicates  that  there  is  no  indication  of  proton 
conduction,  at  least  on  the  time  scale  of  the  simulations.  The 
time  scale  of  "^1  ps  has  been  observed  previously  in  super¬ 
protonic  conductors  as  a  typical  time  scale  associated  with 
the  hydrogen  hopping. The  time  averaged  structure  of 
LiOH  •  2H2O  at  500  K  is  consistent  with  the  (initial)  static 
structure  (Table  II):  all  atoms  including  the  hydrogens  are 
within  ~0.04  A.  The  structure  of  LiOH  •  2H2O  appears  to  be 
stable  at  500  K:  ~50%  of  the  static  internal  coordinates  are 
within  one  standard  deviation  of  each  other  and  the  free  en¬ 
thalpies  of  the  ab  initio  and  the  MD-calculated  phase  agree  to 
within  ~0.6  meV/atom. 

The  absence  of  proton  diffusion  in  these  simulations  in¬ 
dicates  that  the  02-H-03  distance  vectors  do  not  correspond 
to  easy  diffusion  paths  even  at  elevated  temperature.  Thus, 
an  itinerant  network  of  proton  paths  throughout  the  lattice 
does  not  form.  Moreover,  the  simulations  indicate  that  the 
protons  attached  to  the  01  and  02  atoms  do  not  translate. 
This  negative  result  can  be  tentatively  assigned  to  the  poorly 
screened  nuclear  charge  of  the  Li  ions:  the  oxygen  ions  are 
more  strongly  bound,  and  consequently,  the  potential  wells 
for  the  hydrogen  are  deeper  which  makes  proton  hopping  a 
less  likely  event.  The  rather  short  Li-02  distance  of  1.98  A  is 
consistent  with  such  a  strong  influence  of  the  ionic  charge  of 
Li  and  subsequent  confinement  of  the  protons  attached  to  02. 


V.  CONCLUSIONS 

We  report  a  new  dihydrate  of  LiOH  to  form  at  elevated 
pressure.  The  structure  of  this  material  was  determined  by 
comparing  the  ab  initio  calculated  structure  with  the  findings 
from  experimental  synchrotron  x-ray  diffraction  patterns, 
where  calculated  and  observed  structure  factors  agree  with 
an  Rint  of  0.17.  The  structure  exhibits  an  unusual  topology  of 
chains  of  H-bonds  along  (100).  Together  with  a  second  set  of 
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0-H. . .  O  interatomic  distances,  oxygen  and  protons  establish 
an  array  of  Li-filled  channels  along  (101)  compatible  with  the 
Pauling  rules  of  an  ice  network  However,  the  0-H  distances 
of  this  second  set  are  too  long  and  the  O-H4  tetrahedra  too 
distorted  to  establish  H-bonds.  We  examined  the  possibil¬ 
ity  that  elevated  temperature  may  restore  proton  diffusion 
along  those  distance  vectors.  This  would  make  LiOH  •  2H2O  a 
strong  and  lightweight  candidate  for  proton  conduction.  How¬ 
ever,  our  MD  simulations  do  not  show  selective  hydrogen  mo¬ 
tion.  That  is  either  ah  atoms  show  excursions  about  their  equi¬ 
librium  positions  or  ah  atoms  move  simultaneously,  indicat¬ 
ing  that  the  material  is  in  a  liquid  state.  The  disruption  of  the 
H-network  and  the  deep  potential  energy  well  for  the  protons 
both  by  the  poorly  screened  charge  field  of  Li  are  likely  rea¬ 
sons  for  the  absence  of  significant  proton  conduction  in  this 
material. 
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